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ABSTRACT 



The pressure amplitude distribution along the interface between a 
wedge-shaped fluid layer overlying fast and slow fluid bottoms were 
investigated using computer models based on the method of images. 
The pressure amplitude along the bottom rapidly falls off for distances 
closer to the apex than that at which adiabatic-mode theory predicts 
cut-off for the lowest mode. This distance is termed the “dump 
distance” for an acoustically fast-bottom condition. In the case of a 
slow bottom, because the speed of sound ratio is greater than 1 , the 
dump distance does not exist. To facilitate scaling the acoustic fields, 
a “scaling distance” is introduced: Xq = 4 sin 6 c t aiip where 

9c = cos'^ k .2 = (3 is the angle of the wedge, ci = speed of sound 

^2 ^2 

in water, C 2 = speed of sound in bottom, and X = wavelength in the 
wedge medium. Xq appears to have physical meaning as a useful 
scaling distance. 
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I. INTRODUCTION 



Acoustic propagation in an ocean with a sloping bottom is of 
growing concern both for theoretical reasons and because of the 
importance of knowing the performance of acoustic sensors located 
over the continental slope. This problem has been studied by many 
scientists [Ref. 1-81, but many have examined only cases of a fluid 
wedge for which both surfaces are acoustically impenetrable. 

In the case of a fluid wedge for which both surfaces are acousti- 
cally impenetrable, mathematical examination of the problem has been 
performed by Bradley and Hudimic [Ref. 1] using both image theory 
and normal-mode theory. 

For practical interest for underwater acoustics, the top boundary 
may still be approximated as an impenetrable, pressure-release sur- 
face. but the bottom surface can no longer be considered to be per- 
fectly reflecting. From a mathematic point of view, this generalization 
of the boundary conditions has an important repercussion: The wave 
equation is no longer separable and standard normal-mode theory can 
no longer be applied. 

This research will focus on four tasks: 

1. Use the models to reevaluate the fast-bottom case which was 
done by Lesesne [Ref. 9]. 

2. Provide output in graphic form for ease of interpretation. 

3. Change from the fast bottom model to a slow bottom model and 
evaluate whether a scaling factor is useful or not. 

4. Compare results to others’ research. 
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n. DEVELOPMENT 



A. DESCRIPTION OF THE THEORY 
1. Image Theory 

The approach used in this research is to determine the 
amplitude and phase of the pressure in the fluid by the method of 
images (Refs. 8-10). In 1978, with the aid of a computer model. 
Coppens, Sanders, lonnou, and Kawamura [Ref, 11] predicted and 
measured the pressure amplitude and phase in the up-slope direction 
along the bottom of a wedge-shaped fluid layer overlying a fast fluid 
bottom. The geometry is shown in Figure 1. 

Baek [Ref. 12] used the same model in 1984 to predict pres- 
sure amplitude and phase everywhere within the wedge in the up- 
slope direction. In the same year, Lesesne [Ref. 9] implemented a 
model developed by Coppens and Sanders which is not limited to up- 
or down -slope prediction. 

For the wedge-shaped duct, the images which lie on a circle 
(see Figure 2) are eventually going to appear in the region where the 
source is located and which, by the established boundary conditions, 
can contain only the source. However, following Sommerfeld [Ref. 13], 
the concept of an extended Rieman surface is introduced. The exten- 
sion of the 9 coordinate by introduction of Rieman sheets allows for an 
infinite (if need be) number of images. 

The source and each of the images radiate a spherical wave of 
the appropriate phase. The phase coherent summation of these waves 
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yields the total pressure and phase to be found at any field point in the 
wedge. 



The total pressure of all the images and the source along a 



line of constant x from the apex can be determined by the method of 
images (see Figure 2). 



Pn(x) = X (-1)1NT(n/ 2) [gj ^_2 exp[jkix cosOn - D)]] + 



where 

N = lNT(180/p). 

P = the wedge angle. 

D = the receiver angle, 
G = the source angle. 



N 



n=l 



gn exp[jkix cosOn + D)] 



( 1 ) 



= (n - l)p + G, n = 1,3, 5. 7... 



6n = nP - G, 



n = 2, 4, 6,8... 





m= 1,3,5... 





m=2,4,6... 



g -1 = go = +1-0 
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The gn-2 and gn are patJi parameters from the nth image and R(0n) is 
the reflection coefficient of the bottom with grazing angle On- (For 
details, see Ref. 12.) 

The successive reflection of sound from the surface and bot- 
tom is accounted for by multiplying the sound field radiated by each 
image by the plane wave reflection coefficients corresponding to reflec- 
tions encountered by the wave as it propagates toward the field point. 

For the fast bottom, all distances are normalized to the 
distance Xq, measured from the apex, at which the lowest mode would 
attain cut-off in the adiabatic approximation. This distance Xq is 
called the dump distance [Ref. 12] (see Figure 2). 

y ^ 

® tanp 4 * sinOc * tan(i 
X = wavelength in wedge medium 

Q 

0c = cos"^ “ = critical angle 

(3 = wedge angle 

h = water depth at dump distance 

For this model, we use wedge angles of the form n/n where n 
is an integer and consequently no diffraction term exists [Refs. 14-15). 
If n is not an integer, then there will be a diffracted wave which will 
always arrive at a later time than the latest image contribution 
[Ref. 14). 
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2. Normal Mode With Rav Bundles 



V. K. Kuzentsov’s [Ref. 7] normal-mode-with-ray-bundles 
approach treats the entire domain of the wedge as two zones with 
respect to the separation at Vm, where Vm is the exit line of the nor- 
mal mode from the wedge into the half-space. Its position corre- 
sponds to the critical wedge thickness for the given normal mode 
with the non-wave guide zone r < rm. In the wave-guide zone, the 
normal modes behave as normal modes in an ideal wedge. 

3. Rav Tracing 

Tien et al. [Ref. 16) developed an analysis based on a ray 
approach. At the depth where the mode would be cut off, it is 
assumed that a ray reflects from the bottom of the wedge at the criti- 
cal angle: 

This ray is allowed to reflect from the surface and it reattains the 
bottom at a grazing angle which is the critical increased by twice the 
wedge angle. The intervening region of the wedge is filled with a 
large number of downward-traveling rays whose angles of incidence 
on the bottom are uniformly increasing between 0 and 0 + 2|3. Then, 
for each ray the transmission into the bottom is calculated. This 
process is repeated progressively toward the apex of the wedge. 
There is no attempt to preserve phase information, and the power 
of the beam in the bottom at each angle of depression is simply the 
sum of the powers transmitted by the appropriate rays. Calculations 
performed with the help of a computer again provide good agree- 
ment with their experimental observations as far as the angle of 
depression and beam width of the beam in the bottom were con- 
cerned, but this approach predicts that there should be no trans- 
mission of sound into the bottom until the cutoff depth is reached. 
[Ref. 17] 

In the ray-tracing method, there is no attempt to preserve 
phase information, and the power of the beam in the bottom is simply 
the sum of the power transmitted by the appropriate rays. 
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4. Parabolic Equation 

Jensen and Kuperman [Ref. 4] used the parabolic equation 
approximation to study up-slope propagation. Their numerical calcu- 
lations suggested that there is negligible coupling between normal 
modes. They also obtained some very interesting graphs which show 
that the energy of a mode passing through cut-off was radiated as a 
beam into the bottom. 

The only significant limitation of this method appears to be 
that inherent in the approximations necessary for applicability of the 
parabolic equation itself. 

5. Buckingham’s Approach 

In an A. B. Wood Memorial Lecture, Buckingham [Ref. 2] dis- 
cussed acoustic propagation in a wedge-shaped ocean. Basically, he 
follows his previous papers [Ref. 18], modified to include a penetrable 
bottom. In the original paper [Ref. 18], both forward scattering and 
backward scattering exists in the acoustic field except when the 
receiver is located in the vicinity of the source. However, in the 
penetrable wedge there can be no significant backward-reflected 
component in the acoustic field [Ref. 2]. 

Buckingham also used the Rayleigh law of reflection for fast 
bottom and derived an effective wedge with a pressure-release bottom 
(see Figure 31): The Rayleigh law of reflection is 

_ g sin a -Jfcos^a - cos^Oc)^/^ 

“ g sin a + jfcos^a - cos^Oc)^^’^ 
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where g = — , a is the angle from the normal, and a„ is the critical 
Pl ^ 

angle. Here, we are concerned only with grazing angles less than the 

critical, where total internal reflection occurs. When a > ac, there is a 

phase change on reflection, 2E, which from equation (1) can be 

expressed as 



2E = 2 tan~l 



(cos^g - cos^ttc)^/^ 
gsing 



( 2 ) 



Buckingham then points out that similar phase change would occur if 
the ray had been reflected from a pressure-release surface a distance 
A below the actual interface (Figure 31): 




^ksintty 



(3) 



where k = wave number in the first fluid. At the critical angle, it is 
clear from equation (2) that E = 0, in which case. 

A - ^ 

® 2ksinac 



B. DEVELOPMENT OF THE WEDGE MODEL 

Originally, the fast-bottom model was modified from a computer 
model by Coppens, which was programmed in FORTRAN by Baek and 
designed to run on the IBM 3033. The slow-bottom model resulted 
from a modification of this program. For the slow bottom, the cut-off 
condition does not exist. However, using the following formula. 
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Cl = velocity sound in water 



C 2 = velocity sound in bottom 

a scaling factor is introduced which uses the same idea that was in the 
fast-bottom model but with the speed of sound ratio inverted. 

C. DEVELOPMENT OF THE GRAPHIC PROGRAMS 

At this stage, it became necessary to consider alternate forms of 
data presentation. Because of the large volume of computed values, 
numerical output had become ineffective for recognizing details. The 
graphic package Disspla [Ref. 20], available on the IBM 3033 com- 
puter, was used in conjunction with these programs. 

Graphical plots of the data are essential for a clear understanding 
of the problem. A purpose of this research was to utilize graphics to 
better determine the physical concepts necessary to interpret the 
aspects of sound propagation in a wedge-shaped ocean with a pene- 
trable bottom. 

Disspla is an option which uses “function calls” from a FORTRAN 
program to generate graphs. Disspla’s use of FORTRAN allows func- 
tion calls to be written into a program or to be written as separate 
programs. The data used to compute the graphs can, therefore, be 
generated within the FORTRAN program or read in from a data file. 
The latter method was employed since it offered more versatility. 

The size of the data file array or the number of points used did not 
prove to be a limitation on the graphs that were generated. By holding 
all but two parameters constant, a plot of pressure amplitude versus 
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the remaining variable was obtained. The position of the receiver on 
the bottom was usually varied to map the pressure field for a given 
source location. 

The Disspla package supports this, but several problems had to be 
overcome. The Pocket Guide of Disspla describing this option was 
incomplete and very confusing. After comparing the various options 
available, it was decided to use the mode of graphics in which a matrix 
is used to define the data. 

At the beginning, it was necessary to determine the optimal num- 
ber of matrix points which would provide smooth and reasonable plots. 
After running through hundreds of plots, it was decided to use 100 x 
100 matrix points. 

The other important limitation was that the array containing the 
data points must be initially dimensionalized to the exact size of the 
data file array. This information was normally passed from the data file 
during the programs. This made it necessary to initialize the size of 
the data file array before the graphics program was run. 

1. Pressure-Amplitude Graphics 

a The pressure amplitude was plotted as a function of the two 
receiver coordinates. This produced a 3-D plot. The clearest 
presentation was a series of 3-D graphs stacked one after another 
(see Figures 3-A and 4-A), which contain fast- and slow-bottom 
complementary cases. (In this context, “complementary case" 
means the two cases obtained by exchanging the relative values of 
Cl and C 2 .) 

b. When the receiver is close to the source, the spatial variation of 
the pressure amplitude is smaller than the increment, producing 
artificial structure to the graph. To allow a truer view of this 
portion of the pressure amplitude, smaller portions of the field 
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can be viewed with higher resolution without increasing com- 
puter time (see Figures 5 and 6). 

c. When a fuller view was necessary, the models were modified to 
ignore the field close to the source (see Figure 7). 

2. Pressure Contours 

The second graphic program used was a contour plot. The 
following should be mentioned here: 

a The poly 3 interpolation method was used to produce the plots. 

b. The number of matrix points in this program is the same as as for 
the pressure-amplitude graphics. 

c. The increment in the pressure amplitude between contour lines 
was controlled either automatically or by the user. (In Figures 
3-B and 15-B, one can see the difference between the two 
methods.) 

In an attempt to analyze the pressure distribution near the 
source, a polar graphic program was developed. Here, the receiver 
was positioned on circles centered on the source and moved out in 
radial increments. Although this polar form satisfactorily dispersed 
the pressure field near the source, it created very serious geometrical 
distortion (see Figure 8). As a result, it was not very useful. 
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III. RESULTS FROM THE PLOTS 



A. DATA PRESENTATION 

1. The 3-D Plots 

Generally, the source is fixed and the pressure amplitude is 
calculated as a function of receiver position. So, for the three axes of 
the graphs: 

1. The R 2 axis is the normalized receiver distance from the 
shoreline. 

2. The Yq axis is the normalized distance along the shoreline 
between the receiver and the source. 

3. The vertical axis is the pressure amplitude. 

The resolution of the plots is set by specifying a 100 x 100 
matrix of pressure calculation, which means: 

The range of R 2 (or Yp) 

“ The increment of R 2 (or Yq) 

The smaller the range of R2 (or Yq), the smaller the increment of R2 
(or Yq). the higher the resolution of the plots, and the smoother the 
plots. 

2. The Contour Plots 

The axes of the contour plot are X and Y, in which X is equal 
to R 2 and Y is equal to Yq. The label values on the contour lines are the 
pressure amplitudes. The increments of the pressure amplitude and 
the resolution are also controlled by the value of CONMAK [Ref. 20] 
(see Figure 9-B). 
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There are many “tear drop” contours, as shown in Figure 4-B. 
The contours are caused by the limited resolution and disappear at 
higher resolution (see Figure 9-B). 

B. GENERAL PROPERTIES OF THE PLOTS WHICH ARE THE 
SAME FOR BOTH MODELS 

When Ri is equal to an odd number, the pressure peak is directly 
under the source, as shown in Figures 10 and 13. When Ri is equal to 
an even number, the pressure is a minimum under the source (see 
Figures 11 and 14). The explanation for this can be found by using 
only the source, the image character, and the first surface reflected 
path (see Figure 12). 

The time dependence of the signal at a distance di from the 
source is 

Pi = Ae-j(cot-kdi) 



and the time dependence of the signal arriving at R from the image 
characterizing the first surface reflected path is 



p* = ^e-J(wt-kd2) 
1 d2 



where d 2 is the distance from the image point to receiver. So, 



Pt = Pi + Pj = 



Ae-jkdi - Ae-jkd2 
Ldi d2 



= e-J®t 



A_ A g_jk(d2-di) 

di d2 
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Thus, the phase difference between the signals arriving at R from S 
and S' is 

k(d2 - di) = A<1) 

When constructive interference occurs, 

A<J) = 2 tik , n = 0,1,2,... 

and when destructive interference occurs. 

A<|) = (2n + 1 )ti:, n = 0,1,2,... 

The set of points for which d 2 - di has a constant value defines a locus 
along which A({) is constant. The locus d 2 - di = constant is a hyper- 
bola of revolution about the vertical axis connecting the source and the 
first surface-reflected image. The intersection of this surface and the 
plane representing the interface between the two media yields ellipti- 
cal curves close to the source, which evolve into hyperbolic curves far 
from the source. 

These features are observed in the pressure contour plots for the 
various cases. 

The angle of intromission for which the reflection coefficient is 
zero will exist in the fast-bottom case for pi > p 2 or the slow-bottom 
case for pi < p 2 . Figures 17 and 18 are for cases possessing angles of 
intromission; these appear to exhibit no significant pressure changes 
compared to the non-intromission conditions. 
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C. DIFFERENCE BETWEEN BOTH MODELS 

With the fast bottom, the pressure pattern is much more irregular 
and the pressure decays with distance from the source less rapidly 
because, with the fast bottom, more images are totally reflected from 
the bottom. 
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IV. EVALUATION 



The effects of changing the properties of the wedge and the 
source position on the prediction of the pressure field were evaluated. 
As a base-line case, the following parameters were used: 
p = 10° 

G =5° 

D =0° 

Rl = 40 

^ = 0.9 

P2 

Ci/C 2 = 0.9 (fast bottom) 

ci/C 2 = 1/0.9 (slow bottom) 

a/k2 = 0.0001 

Three-dimensional depictions of the pressure fields for these 
conditions as functions of R 2 and Yq are shown in Figures 3 and 4 for 
both fast- and slow-bottom cases. The corresponding contour plots 
are included. 

Both models were tested by varying one parameter and holding 
the others constant. Comparison of pressure-contour plots as well as 
three-dimensional plots was then used to determine the effect of each 
parameter. 

A. SPEED-OF-SOUND RATIO, DENSITY RATIO, AND ATTENUATION 

The base-line case was first examined by changing the speed-of- 
sound ratio. The positions and number of crests exhibited simple and 
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clear patterns for the slow-bottom model (see Figure 19) and more 
complex patterns for the fast-bottom model (see Figure 20) as this 
ratio was changed from 0.9 (or 1/0.9) to 0.5 (1/0.5). The density ratio 
was then changed down to a value of 0.5. The positions and number of 
crests show only minor variation (see Figure 21). This suggests that 
the major features of the sound field in the wedge are relatively insen- 
sitive to minor variations in the characteristics of the bottom. 

Finally, the value of a/k 2 was increased from 0.0001 to 0.1. The 
pressure-pattern did not change significantly but, as expected, the 
amplitudes of the crests showed sharp decreases as the attenuation 
increases. 

B. RECEIVER ANGLE 

The reader should be aware of the fact that the program con- 
tained in the appendix was modified as simply as possible, to insure 
that pressure was calculated infinitesimally above the bottom. If it is 
desired to calculate the pressure for receiver angles not essentially 
zero, the leading statements and do loops must be modified. 

C. WEDGE ANGLE 

As the wedge angle decreases, the dump distance increases, the 
number of images increases, and therefore the CPU time increases. 
For example, for P = 10°, the CPU times for fast-bottom and slow-bot- 
tom cases are about 520 seconds and 480 seconds, respectively. If P = 
7°, the CPU time increases to about 900 seconds (see Figures 22 and 
23). 
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D. SOURCE ANGLE 

Changing the distance of the source from the apex (but still with 
G = p/2) had a great effect on the position of the crests and troughs. 
However, for the source at large distances from the apex, scaling the 
receiver position restores the original patterns. For example, if the 
source is moved outward to twice its original position, a plot with the 
length scales divided by two produces the same graphic presentation 
(see Figures 3 and 24 or 4 and 25). This is consistent with Lesesne’s 
observation. 

E. LLOYD’S MIRROR 

Figures 29 and 30 show the field due to the source and the first 
image point only. The source and the first image are of equal strength. 
The first image is coherent with the source but k out of phase (see 
Figure 12). That produces the Lloyd’s mirror. The pressure 
distribution observed is identical for both fast- and slow-bottom 
programs when they were modified to have the same apparent speed 
of sound in the wedge. These figures show that major features of the 
pressure field are due mainly to the source and first image when the 
receiver is close to the source. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The theory of 3-D acoustic propagation in a penetrable wedge 
presented in this research is based on the method of images which 
includes the fast- and the slow-bottom models. One of the research 
goals was to test whether the scaling factor method can apply to the 
slow-bottom case. 

A. CONCLUSIONS 

For both the fast bottom and the slow bottom, when the receiver 
is near the source, strong spatial oscillations in pressure are present. 
These variations are mainly due to interference between the source 
and the first surface-reflected image. 

The primary result of the fast-bottom model is the importance of 
interference evident in the form of strong spatial oscillation. When 
the bottom mismatch increases, more interference occurs, because 
more images contribute to the field. 

Generally, investigation of this model showed that the scaling fac- 
tor worked. 

One observation is in violation of physical intuition— when the 
bottom mismatch increases, more interference should occur. Figures 
26. 27, and 28 show the opposite, with less interference as the bot- 
tom mismatch increases. 

The results for the fast-bottom model are more complex than for 
the slow-bottom case because the cut-off condition only exists in the 
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fast-bottom model, which yields more images which correspond to 
perfect reflections for the bottom. 

The method of images does not yield an easy decomposition of the 
sound field into a form that can be closely associated with normal 
modes. On the other hand, the method of images provides intrinsi- 
cally the phase interference effects associated with the superposition 
of effects that would be more difficult to extract from the adiabatic 
normal mode. The intrinsic inclusion of absorption in the substrate is 
an attractive and realistic feature of the image approach. In addition, 
the image approach appears to avoid the generation of caustics at 
turning points which are inextricably linked with the results of 
unmodified ray tracing approaches. 

B. RECOMMENDATIONS FOR FURTHER INVESTIGATION 

1. For the case of a real fast-bottom, identification of the normal 
modes of the system and investigation of the possibility of normal 
mode coupling is worth study. 

2. Using the batch system to investigate the small wedge angle is 
worth study. 

3. Changing the position and depth of the receiver is worth study. 

4. For the slow bottom, experiments should be undertaken to com- 
pare with the theoretical predictions. 
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APPENDIX A 



COMPUTER PROGRAMS 



FILE: 3LS4 FORTRAN A1 



DEBUG UNIT(6),SUBCHK 
END DEBUG 



CROSS SLOP TRANSMISSION IN THE WEDGE WITH PENETRABLE LOSSY SLOW 
BOTTOM ( NO GEOMETRICAL APPROXIMATIONS ) 



DIMENSIONZ(M.M) (21,11), (Y0,R2)--(L.A) 

DIMENSION 2(101,101) 

COMMON WORK(SOOOO) 

INTEGER A,II,JJ,N1,S1,S2,IE,N,I 

REALM RR1 ,R2,Y0,Y1 ,YOSTA,YOINC,R2INC,R2STA,RZ,RZ1 ,M,J, 

1 B,C(30).CC,C2,D,D1 ,D2,E(30),E1 ,E2,F(30),F1 ,F2,G,H,PI,FI, 

1 PI ,P2,P(1 25,1 30),PH,Q1 ,Y2,R3,R8(80),R9(80), 

1 S(30),T,T4,T6,T1 (80),TQQ,TQQl ,TQQ2,TQQ3,W,W0, W1 ,XL,Y, 

1 Z1,Z2,Z3,Z4,Z5,Z6,Z,EL,DIST 

CALL SETIME 
PI = ARCOS(-I.O) 

c ..................... 

C INPUT PARAMETERS 

C N1 = # OF IMAGINE POINTS 

C B = WEDGE ANGLE (DEG) 

C G = SOURCE ANGLE (DEG) 

C D = RECEIVER ANGLE (DEG) 

C R1 = SOURCE DISTANCE (IN DUMP DISTANCES) 

C R2 = RECEIVER DISTANCE (IN DUMP DISTANCES) 

C YO = SHORE DISPLACEMENT (IN DUMP DISTANCES) 

C D1 = RH01/RH02 

C CC = C1/C2 

C XL = ALPHA/K2 

C A = # OF RECEIVER POSITIONS IN R2 DIRECTION 

C L = # OF RECEIVER POSITIONS IN YO DIRECTION 

C R2INC = INCREMENTAL INCREASE FOR R2 

C Y0INC = INCREMENTAL INCREASE FOR YO 

C 

Q ********************************* ******0y^y‘y^ ^^1^0 ^2** ««*«««*«***«********«* ********************** 

C CALL TEK618 

CALL COMPRS 
C CALL PAGE(1 1,8.5) 

B=10. 

G = 1.0*B/4. 
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c 

c 



G = 5 
D = 0.0 
R1 = 10. 

R2 = 5. 

Y0 = 0. 

D1 = .9 
D1 = 1/.9 

FOR SLOW-BOTTOM CASE 
CC = 1/0.9 
C FOR FAST-BOTTOM CASE 

CC = 0.9 
C CC = .5 

XL = .0001 
A = 100 
L= 100 
R2INC = .1 
Y0INC = .1 

Q *»***•**••*##•*•#******** *********«*^y^j^ *«*•** ****«***•*«**««**•**« ****•« 

YOSTA = YO 
R2STA = R2 
W = 0.0 

WRITE(6,300) B.Q.D 

300 FORMAT (3X.’B = •.F5.2,3X,’G = '.F5.2,3X,’D = ’.F5.2.3X) 

310 FORMAT (1X.F5.2,1X,F5.2,1X,F5.2) 

WRITE(6,400) R1 

400 FORMAT (3X, ’SOURCE RANGE = ’,F8.2) 

410 FORMAT (1X.F6.2) 

WRITE(6,500) Dl.CC.XL 

500 FORMAT (3X,’DENSITY RATIO = ’.F8.6.3X, ’SPEED RATIO = ’,F8.6,3X, 

* ’ALPHA/K2 = ’,F8,6) 

510 FORMAT (1X,F6.4,1X,F6.4,1X,F6.4) 

WRITE(6,600) 

600 FORMAT (5X,’42’,6X,’Y0’,1 OX.’P') 

C WRITE(22,700) R2.Y0 
700 FORMAT (IX, F6.2, IX, F6.2) 

N1 = (180/B + .00001) 

T6 = 180/PI 
B = B/T6 
G = G7T6 
C2 = CC**2 
C 

HM=0 

HH=1 

D =HM*B/A 
DD = B/A 
V = 2*B/10 

IF(D.LT.V) GO TO 110 
IF(D.GE.V) GOTO 120 
110 DX = D 

C THIS DO LOOP CALCULATES THE TMETA(N) AND THE IMINAGE SL4NT 
C RANGES R8(N) AND R9(N). 

C 
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c 



DO 55 HH = 1.10 
D = DX+(HH-1)*B/(10*A) 

DD = B/(10*A) 

120 S1=1 

T4 = P1/(2*TAN(ARCOS(1/CC))*TAN(B)) 
TOO = TAN(B) 

FOR SLOW-BOTTOM CASE 
TOOi = ARCOS(1/CC) 

FOR FAST-BOTTOM CASE 
TOOi = ARCOS(CC) 

TOQ2 = TAN(TOOI) 

TOQ3 = 2.*TQ02*T0Q 
T4 = PI/TQQ3 
01 = 1/SORT(2.0) 

LP1X=L+1 
IAP1X=A+1 
DO 20 M=1,LP1X 
DO10J=1,IAP1X 
D2 = (Y0**2)+(R1**2)+(R2**2) 

R3 = 2*R1*R2 



DO 30 N = 1,N1 

IF(SI.GT.O) T1(N)=FLOAT(N-1)‘B+G 
IF(SI.LT.O) T1(N)=FLOAT(N)*B-G 
SI =-S1 

R8(N) = SQRT(D2-R3*COS(T1(N)-D)) 
R9(N) = SORT(D2-R3*COS(T1(N)+D)) 
CONTINUE 



SUM PRESSURE OVER ALL IMINAGES- 



PI = 0.0 
P2 = 0.0 
DO 40 N = 1,N1 

S2 = (-1)**(IFIX(N/2+.0001)) 
W1 = 2*C2*XL 



FEFL COEFFS ALONG WITH UPPER PATH— 



II = IFIX((N-.9999)/2) 

DO 41 1 = 1,11 

S(l) = ABS(R1*SIN(T1(N)-2*FLOAT(l)*B) 
+R2*SIN(2*FL0AT(I)*B-D))/R8(N) 
C(l) =SQRT(1-(S(I)**2)) 

T = S(l)/D1 

WO = (-C2+(C(I)**2)) 

Y = SQRT((W0**2)+(W1**2)) 

Z(M,J) = ABS(WO) 

IF(Y,LE,Z(M,J)) Y = Z(M,J) 

Y1 = Q1*SQRT(Y+W0) 
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OOO OOOO OOO 0^000 



Y2 = -Q1*SQRT(Y-W0) 

E(l) = (T**2)-(Y2**2)-((Y1**2))/(((T-Y2)**2)+(Y1**2) 
F(l) = (2*Y1*T)/(((T-Y2)**2)+(Y1**2)) 

41 CONTINUE 



PRODUCT OF REFL COEFFS ALONG NTH UPPER PATH 



El = 1.0 
FI = 0.0 

IF(N.LE.2.0) GOTO 422 
DO 42 1 = 1,11 

E2 = E1*E(I)-F1*F(I) 
F2 = F1*E(I)-E1*F(I) 
El = E2 
FI = F2 
CONTINUE 



FOR SLOW-BOTTOM CASE 
T = T4*R8(N) 

FOR FAST-BOTTOM CASE 
T = T4*Rs(N)/CC 



PI =P1 +S2*(E1*COS(T)+F1*SIN(T))/48(N) 
P2 = P2 + S2*(F1*C0S(T)+E1*SIN(T))/48(N) 



REFL COEFFS ALONG NTH LOWER PATH 



II = II +1 
D043 1 = 1,11 

S(l) = ABS(R1*SIN(T1(N)-2*FLOAT(l-1)‘8) 

+ R2*SIN(2*FLOAT(M )*B+D))/R9(N) 

C(I) = SQRT(1-(S(I)**2)) 

T = S(I)/D1 

WO = (-C2+(C(I)*‘2)) 

Y = SQRT((W0**2)+(W1**2)) 

Z(M,J) = ABS(WO) 

IF(Y.LE.Z{M,J))Y = Z(M,J) 

Y1 = Q1*SQRT(Y+W0) 

Y2 = -Q1*SQRT(Y-W0) 

E(l) = ({T**2)-(Y2**2)-(Y1**2))/(((T-YW)**2)_(Y1**2)) 
F(l) = (2*Y1*T)/(((T-Y2)**2)+(Y1**2)) 

AN2 = TAN(F(I)/E(I))*T6) 

CONTINUE 



PRODUCT OF REFL COEFFS ALONG NTH LOWER PATH 



El = 1. 

FI =0. 

DO 44 1 = 1,11 

E2 = E1*E(I)-F1*F(I) 
F2 = F1*E(I)+E1*F(I) 
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El = E2 
F1 =F2 

44 CONTINUE 

0 ♦*******♦♦♦***♦***♦*♦♦**♦*♦♦♦***♦****♦****♦******♦*********♦*♦*♦ 
C FOR SLOW-BOTTOM CASE 

T = T4*R9(N) 

C FOR FAST-BOTTOM CASE 

T = T4*Rg(N)/CC FAST BOTTOM 
0 ************«**********««*«««*««***«*«*******«««***«*****«****«* 
P1 = P1 + S2*(E1*COS(T)+F1*SIN(T))/R9(N) 
P2 = P2 + S2*(F1*C0S(T)-E1*SIN(T))/R9(N) 
40 CONTINUE 

Z(M,J) = SQRT((P1“2)+D(P2‘*2))*R1 
IF(W.LT.Z(M,J)) W=Z(M,J) 

WRITE(6,200) R2,Y0,Z(M,J) 



200 FORMAT (3X,F7.3,3X,F7.2,3X,F9.5,3X,F9.5) 

C 

0 **♦♦*♦♦*♦♦******♦♦**♦♦♦♦♦**♦*****♦♦♦♦********♦**♦*♦♦*♦*♦♦*♦♦♦****♦♦*♦♦*♦♦ 

R2 = R2+R2INC 

Q 

C DIST = SQRT((J-41)**2) + (M-1)**2) 

C IF (DIST.LE.20.) Z(M,J) = 0. 

C * 

10 CONTINUE 

Y0 = Y0 + Y0INC 
R2 = 5. 

20 CONTINUE 

CALLSURF3D(Z,NX,NY) 

YO = 0. 

55 CONTINUE 

33 CONTINUE 

STOP 
END 

C 

SUBROUTINE SURF3D(Z,NX,NY) 

C 

C PLOTS 3D SURFACE OF AN ARRAY Z(NX,NY) OF REGULARLY SPACED POINTS. 

C 

DIMENSION Z(101,101) 

DIMENSION RZ(IOI.LIOI) 

COMMON WORK(50000) 

C 

C CHECK DATA STATEMENTS TO ET AXIS LIMITS; 

C 

DO 1 1-1,101 
DO 1 J-1,101 
1 RZ(J,I)=Z(I,J) 

c 

C PRODUCE CONTOUR PLOT 

C 

CALL PAGE (11,8.5) 
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CALL SHDCHR(.90,1,, 003,1.) 

CALL AREA2D(9.5.6) 

C CALL NOBRDR 

CALL BLOWUP(.85) 

CALL XNAME(’XXXM) 

CALL YNAME(’YYYM) 

CALL HEADINCCONTOUR PLOT$MOO,1.,1) 
CALL GRAF(5.,1.,15.,0.,1.,10.) 

CALL FRAME 
CALL BCOMON(50000) 

CALL CONMAK(RZ,101,101,.5) 

CALL CONLIN(0.’SOLID’,'LABELS’,1 ,1) 
CALL CONANG(90.) 

CALL POLY3 

CALL CONTUR(2, ’LABELS’, -DRAW’) 

CALL ENDPL(O) 



PRODUCE SURFACE PLOT- 



CALL PAGE(11,8.5) 

CALL AREA2D(9.5,6) 

CALL MX1ALF(’STANDARD’,’*’) 
CALL MX2ALF(’UCSTD’,’-h’) 



LABEL AXES: 

CALL X3NAME(’R2 A-fXIS‘$’,100) 

CALL Y3NAME(’Y0 A+XIS*$’,100) 

CALL Z3NAME(’PRESSURE A-fXIS*$’,100) 



CALL VOLM3D(5.,5.,5.) 

CALL YINTAX 
CALL XINTAX 
CALL ZINTAX 
CALL ZAXANG(0.) 

CALL YAXANG(0.) 

CALL XTICKS(5) 

CALL YTICKS(5) 

CALL HEADIN (’SLOW BOTTOM’,1 1,1.5,1) 

CALL HEADIN (’FAST BOTTOM’,1 1 ,1.5,1) 

CALL MESSAG (’C.P.U. = SECONDS$’,100,6.5,0.) 

CALL MESSAG (’VIEPTS IS 250,350,90$’,1 00.6,6.5) 

CALL MESSAG (’SURMAT IS RZ,1,101,101,101$’,100,6,6) 

CALL MESSAG ('SURVIS IS TOP$’,10,6,5.5) 

CALL MESSAG {’B(W.A.)-10.,G(S.A)-5.,D(R.A)=0 $’,100,6,5) 

CALL MESSAG (’RH01/RH02 =.9,ALPHA/K2=.0001 $’,100,6.4.5) 
CALL MESSAG (’C1/C2 -1/.9,R1 = 10.,42=5.,Y0=0.$’,100,6,4) 
CALL MESSAG (’R2INC =.1,Y0INC =.1,CMK =.5$’,100,6,3.5) 



CALL VUABS(250.,350.,90.) 

CALL SURVIS(’TOP’) 

CALL GRAF3D (5.,1.,15.,0.,1.,10.,0.,5.,15.) 
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ALL BLSUR 

C CALLSURMAT(RZ,XX,A+1,XX,L+1) 

CALL SURMAT(RZ, 1,101, 101, 101,0) 

C 

CALL ENDPL(O) 

CALL DONEPL 
CALL GETIME(IET) 

EL = lET * .000026 
WRITE (4,2) EL 

2 FORMAT(20X,TIME = ’,E16.6.’SECONDS.') 

RETURN 
END 
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FIGURES 
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Figure 1 

Three-Dimensional Wedge Geometry 




Figure 2 

The Geometry of Image Solution 



40 



'T-r 







< 

I 

CO 

o 



'ZXj 



41 



Pressure on the Bottom of the Wedge as Calculated by the Image Model 
for the Indicated Values of Parameters 




42 



Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 



SURriRT IS RZ, 1,101,101,101 
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Pressure on the Bottom of the Wedge as Calculated by the Image Model 
for the Indicated Values of Parameters 
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Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 



suRimi 
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Pressure on the Bottom of the Wedge as Calculated by the Image Model 
for the Indicated Values of Parameters 
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Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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Pressure on the Bottom of the Wedge as Calculated by the Image Model 

for the Indicated Values of Parameters 
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Figure 6-B 

Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 









o 




o 










ii 




a 


• 








- — . 


» — < 


o 


» — I 








cz 


O 


>- 










• 


o 




B 








CO 


o 


o 






O 




< — ^ 


• 


8 






t— 1 




c 


B 














CM 


CM 


CJ 




, — 1 




• 


CO 


CO 






o 




LO 






• 


o 








cr 


o 


t — 1 


CD 






a 


~JZ 






•s 


T— « 




^ — • 


Q_ 




B 


o 


o 




ci: 


_I 


fl 




LO 






• 


cn 




O 








CO 




— 1 


~T^ 


•v 


, — 1 




— ' 


CD 


CO 


• » 


o 






CO 


• 




o 


LO 




o 


-- 


O 


cn 


>- 


OJ 


a: 


t — ' 


• 


H 


• 










o 






• 


oo 


00 


CO 


, — . 


OO 


— 1 


^ — I 


• — 1 


1 — . 


. — . 


a 


o 

















II 


a 


CO 


h- 


CO 


• 


CO 






e— 


cn 


► — 1 


n~ 




oj 


CJ 


n 


20 




• 


, — < 


o 


~z 


LJ 


DO 


CO 




o 




}—4 


. — I 


Z) 


r: 




JZ 


, — . 


CM 


> 


CO 


CO 


CQ 


CO 


CJ 


CO 






1 ^ 

01 ^ 



CO 

g 

o 

CJ 

LiJ 

Lll 




Q 



49 



Figure 7-A 

Pressure on the Bottom of the Wedge as Calculated by the Image Model 

for the Indicated Values of Parameters 




50 



Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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Pressure on the Bottom of the Wedge as Calculated by the Image Model 

for the Indicated Values of Parameters 
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Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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Pressure on the Bottom of the Wedge as Calculated by the Image Model 

for the Indicated Values of Parameters 
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Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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Pressure on the Bottom of the Wedge as Calculated by the Image Model 

for the Indicated Values of Parameters 
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Figure 10-B 

Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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Pressure on the Bottom of the Wedge as Calculated by the Image Model 

for the Indicated Values of Parameters 
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Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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S = Source point 
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s' = Image point 



Figure 12 

Interference Pattern from a Point Source and an Image Point 
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Figure 13 A 

Pressure on the Bottom of the Wedge as Calculated by the Image Model 
for the Indicated Values of Parameters 
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Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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Pressure on Axis Response of Wedge 
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Pressure on the Bottom of the Wedge as Calculated by the Image Model 

for the Indicated Values of Parameters 
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Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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Pressure on the Bottom of the Wedge as Calculated by the Image Model 

for the Indicated Values of Parameters 
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Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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Figure 17-A 

Pressure on the Bottom of the Wedge as Calculated by the Image Model 
for the Indicated Values of Parameters 
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Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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Pressure on the Bottom of the Wedge as Calculated by the Image Model 

for the Indicated Values of Parameters 
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Figure 18-B 

Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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Pressure on the Bottom of the Wedge as Calculated by the Image Model 

for the Indicated Values of Parameters 
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Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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Pressure on the Bottom of the Wedge as Calculated by the Image Model 

for the Indicated Values of Parameters 
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Eiqual-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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Pressure on the Bottom of the Wedge as Calculated by the Image Model 

for the Indicated Values of Parameters 
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Figure 21-B 

Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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Figure 22-A 

Pressure on the Bottom of the Wedge as Calculated by the Image Model 

for the Indicated Values of Parameters 
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Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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Pressure on the Bottom of the Wedge as Calculated by the Image Model 

for the Indicated Values of Parameters 
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Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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Pressure on the Bottom of the Wedge as Calculated by the Image Model 

for the Indicated Values of Parameters 
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Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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Figure 25-A 

Pressure on the Bottom of the Wedge as Calculated by the Image Model 

for the Indicated Values of Parameters 
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Figure 25-13 

Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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ccNiouR Plot 
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Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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Pressure on the Bottom of the Wedge as Calculated by the Image Model 
for the Indicated Values of Parameters 
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Pressure on the Bottom of the Wedge as Calculated by the Image Model 
for the Indicated Values of Parameters 
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Figure 29-B 

Equal-Pressure Contours on the Bottom of the Wedge 
for the Indicated Values of Parameters 
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Pressure on the Bottom of the Wedge as Calculated by the Image Model 
for the Indicated Values of Parameters 
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Values of Parameters 





(b) 



(a) Total internal reflection of a ray and 

(b) a reflection with the same phase shift from an 
effective-release surface 




(a) Wedge with penetrable bottom and apex 0 
and "effective” wedge wth pressure-release 
bottom and apex 0. 



Figure 31 

Effective Pressure-Release Bottom 
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